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1 Introduction 

Aspects of gauge theory, Hamiltonian mechanics, relativity and quantum mechanics arise natu- 
rally in the mathematics of a non-commutative framework for calculus and differential geome- 
try. In this paper, we first give a review of our previous results about discrete phyics and non- 
commutative worlds. The simplest discrete system corresponds directly to the square root of 
minus one, seen as an oscillation between one and minus one. This way thinking about i as an 
iterant is explained below. By starting with a discrete time series of positions, one has immedi- 
ately a non-commutativity of observations since the measurement of velocity involves the tick of 
the clock and the measurment of position does not demand the tick of the clock. Commutators 
that arise from discrete observation suggest a non-commutative calculus, and this calculus leads 
to a generalization of standard advanced calculus in terms of a non-commutative world. In a non- 
commutative world, all derivatives are represented by commutators. We then give our version 
of Feynman-Dyson derivation of the formalism of electromagnetic gauge theory. The rest of the 
paper investigates algebraic constraints that bind the commutative and non-commutative worlds. 

Section 2 is a self-contained miniature version of the whole story in this paper, starting with 
the square root of minus one seen as a discrete oscillation, a clock. We proceed from there and 
analyze the position of the square root of minus one in relation to discrete systems and quantum 
mechanics. We end this section by fitting together these observations into the structure of the 
Heisenberg commutator 

[p, q] = ih. 

Section 3 is a review of the context of non-commutative worlds with a preview of part of the 
Feynman-Dyson derivation. This section generalizes the concepts in Section 2 and places them 
in the wider context of non-commutative worlds. The key to this generalization is our method of 
embedding discrete calculus into non-commutative calculus. Section 4 is a discussion of iterants 



and matrix algebra. We show how matrix algebra in any dimension can be regarded as describ- 
ing the pattern of acts of observation (time shifting operators corresponding to permutations) on 
periodic time series. Section 5 gives a complete treatment of our generalization of the Feynman- 
Dyson derivation of Maxwell's equations in a non-commutative framework. This section is the 
first foray into the consequences of constraints. This version of the Feynman-Dyson derivation 
depends entirely on the postulation of a full time derivative in the non-commutative world that 
matches the corresponding formula in ordinary commutative advanced calculus. Section 6 dis- 
cusses constraints on non-commutative worlds that are imposed by asking for correspondences 
between forms of classical differentiation and the derivatives represented by commutators in a 
correpsondent non-commutative world. This discussion of constraints parallels work of Tony 
Deakin (31 and will be continued in joint work of the author and Deakin. At the level of 
the second constraint we encounter issues related to general relativity. Section 7 continues the 
constraints discussion in Section 5, showing how to generalize to higher-order constraints and 
obtains a commutator formula for the third order constraint. The first appendix, Section 8, is a 
very condensed review of the relationship of the Bianchi identity in differential geometry and 
the Einstein equations for general relativity. We then observe that every derivation in a non- 
commutative world comes equipped with its own Bianchi identity. This observation suggests one 
way to investigate general relativity in the non-commutative context. Section 9 is a philosophical 
appendix. 



2 Quantum Mechanics - The Square Root of Minus One is a 
Clock 

The purpose of this section is to place i, the square root of minus one, and its algebra in a 
context of discrete recursive systems. We begin by starting with a simple periodic process that is 
associated directly with the classical attempt to solve for i as a solution to a quadratic equation. 
We take the point of view that solving x 2 = ax + b is the same (when x ^ 0) as solving 

x = a + b/x, 

and hence is a matter of finding a fixed point. In the case of i we have 

x 2 = -1 

and so desire a fixed point 

x = —1/x. 

There are no real numbers that are fixed points for this operator and so we consider the oscillatory 
process generated by 

R(x) = -1/x. 

The fixed point would satisfy 

% = -1/i 
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+1, -1, +1, -1, +1, -1, +1, -1 



[-1.+1] [+1.-1] 

Figure 1 : A Basic Oscillation 

and multiplying, we get that 

ii = — 1. 

On the other hand the iteration of R yields 

1, R(l) = -1, R(R(1)) = +1, R(R(R(1))) = -1, +1,-1, +1, -1, • • • ■ 

The square root of minus one is a perfect example of an eigenform that occurs in a new and wider 
domain than the original context in which its recursive process arose. The process has no fixed 
point in the original domain. 

Looking at the oscillation between +1 and —1, we see that there are naturally two phase- 
shifted viewpoints. We denote these two views of the oscillation by [+1, —1] and[— 1, +1]. These 
viewpoints correspond to whether one regards the oscillation at time zero as starting with +1 or 
with —1. See Figure 1. 

We shall let J{+1, —1} stand for an undisclosed alternation or ambiguity between +1 and 
— 1 and call 7{+l, —1} an iterant. There are two iterant views: [+1, —1] and [—1, +1]. 

Given an iterant [a, b], we can think of [b, a] as the same process with a shift of one time step. 
These two iterant views, seen as points of view of an alternating process, will become the square 
roots of negative unity, i and —i. 

We introduce a temporal shift operator 77 such that 

[a,b]r) = rj[b, a] 

and 

7777 = 1 

for any iterant [a, b], so that concatenated observations can include a time step of one-half period 
of the process 

■ ■ • abababab ■ ■ ■ . 
We combine iterant views term-by-term as in 

[a, b] [c, d] = [ac, bd}. 
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II = 



[1, -l] v = [1, -!][-!, l] w = [-1, -1 



1. 



Thus we have constructed the square root of minus one by using an iterant viewpoint. In this 
view i represents a discrete oscillating temporal process and it is an eigenform for R(x) = — 1/x, 
participating in the algebraic structure of the complex numbers. In fact the corresponding algebra 
structure of linear combinations [a, b] + [c, d\q is isomorphic with 2x2 matrix algebra and iterants 
can be used to construct n x n matrix algebra. We treat this generalization in Section 4 of this 
paper. 

The Temporal Nexus. We take as a matter of principle that the usual real variable tfor time is 
better represented as it so that time is seen to be a process, an observation and a magnitude all at 
once. This principle of "imaginary time" is justified by the eigenform approach to the structure 
of time and the structure of the square root of minus one. 

As an example of the use of the Temporal Nexus, consider the expression x 2 + y 2 + z 2 + t 2 , 
the square of the Euclidean distance of a point (x, y, z, t) from the origin in Euclidean four- 
dimensional space. Now replace t by it, and find 



the squared distance in hyperbolic metric for special relativity. By replacing t by its process 
operator value it we make the transition to the physical mathematics of special relativity. 

2.1 Quantum Physics, Eigenvalue and Eigenform 

In quantum modeling, the state of a physical system is represented by a vector in a Hilbert space. 
As time goes on the vector undergoes a unitary evolution in the Hilbert space. Observable quanti- 
ties correspond to Hermitian operators H and vectors v that have the property that the application 
of H to v results in a new vector that is a multiple of v by a real factor A. Thus 



One says that v is an eigenvector for the operator H, and that A is the eigenvalue. 

2.2 The Wave Function in Quantum Mechanics and The Square Root of 
Minus One 

One can regard a wave function such as ip(x,t) = exp(i(kx — wt)) as containing a micro- 
oscillatory system with the special synchronizations of the iterant view i = [+1, —l]r) . It is these 



x 2 + y 2 + z 2 + {it) 2 



x 2 + y 2 + z 2 - t 2 , 



Hv = Xv. 
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synchronizations that make the big eigenform of the exponential work correctly with respect to 
differentiation, allowing it to create the appearance of rotational behaviour, wave behaviour and 
the semblance of the continuum. In other words, we are suggesting that once can take a temporal 
view of the well-known equation of Euler: 

e %e = cos(6) + isin{6) 

by regarding the i in this equation as an iterant, as discrete oscillation between —1 and +1. One 
can blend the classical geometrical view of the complex numbers with the iterant view by thinking 
of a point that orbits the origin of the complex plane, intersecting the real axis periodically and 
producing, in the real axis, a periodic oscillation in relation to its orbital movement in the two 
dimensional space. The special synchronization is the algebra of the time shift embodied in 

7777 = 1 

and 

[a, % = rj[b, a] 

that makes the algebra of i = [1, — l]r) imply that i 2 = — 1. This interpretation does not change 
the formalism of these complex- valued functions, but it does change one's point of view and we 
now show how the properties of i as a discrete dynamical systerm are found in any such system. 

2.3 Time Series and Discrete Physics 

We have just reformulated the complex numbers and expanded the context of matrix algebra to an 
interpretation of i as an oscillatory process and matrix elements as combined spatial and temporal 
oscillatory processes (in the sense that [a, b] is not affected in its order by a time step, while [a, b]r] 
includes the time dynamic in its interactive capability, and 2x2 matrix algebra is the algebra of 
iterant views [a, b] + [c, djrj). 

We now consider elementary discrete physics in one dimension. Consider a time series of 
positions 

x(t) : t = 0, At, 2At, 3At, 
We can define the velocity v(t) by the formula 

v{t) = (x(t + At) - x(t))/At = Dx{t) 

where D denotes this discrete derivative. In order to obtain v(t) we need at least one tick At of 
the discrete clock. Just as in the iterant algebra, we need a time-shift operator to handle the fact 
that once we have observed v(t), the time has moved up by one tick. 

We adjust the discrete derivative. We shall add an operator J that in this context accomplishes 
the time shift: 

x{t)J = Jx{t + At). 
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We then redefine the derivative to include this shift: 

Dx{t) = J(x(t + At) - x(t))/At. 

This readjustment of the derivative rewrites it so that the temporal properties of successive obser- 
vations are handled automatically. 

Discrete observations do not commute. Let A and B denote quantities that we wish to observe 
in the discrete system. Let AB denote the result of first observing B and then observing A. The 
result of this definition is that a successive observation of the form x(Dx) is distinct from an 
observation of the form (Dx)x. In the first case, we first observe the velocity at time t, and then 
x is measured at t + At. In the second case, we measure x at t and then measure the velocity. 
We measure the difference between these two results by taking a commutator 

[A, B] = AB- BA 

and we get the following computations where we write Ax = x(t + At) — x(t). 

x(Dx) = x(t) J(x(t + At) - x(t)) = Jx{t + At)(x(t + At) - x(t)). 
(Dx)x = J(x(t + At) - x(t))x(t). 
[x, Dx] = x(Dx) - (Dx)x = (J/At)(x(t + At) - x(t)) 2 = J(Ax) 2 /At 
This final result is worth recording: 

[x, Dx] = J(Ax) 2 /At. 

From this result we see that the commutator of x and.Dx will be constant if (Ax) 2 /At = K is a 
constant. For a given time-step, this means that 

(Ax) 2 = KAt 

so that 

Ax = ±y/(KAt) 
This is a Brownian process with diffusion constant equal to K. 

Thus we arrive at the result that any discrete process viewed in this framework of discrete 
observation has the basic commutator 

[x,Dx] = J(Ax) 2 /At, 

generalizing a Brownian process and containing the factor (Ax) 2 /At that corresponds to the 
classical diffusion constant. It is worth noting that the adjusment that we have made to the 
discrete derivative makes it into a commutator as follows: 

Dx(t) = J(x(t + At) - x(t))/At = (x(t)J - Jx(t))At = [x(t), J]/ At. 

By replacing discrete derivatives by commutators we can express discrete physics in many vari- 
ables in a context of non-commutative algebra. We enter this generalization in the next section 
of the paper. 

We now use the temporal nexus (the square root of minus one as a clock) and rewrite these 
commutators to match quantum mechanics. 
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2.4 Simplicity and the Heisenberg Commutator 



Finally, we arrive at the simplest place. Time and the square root of minus one are inseparable 
in the temporal nexus. The square root of minus one is a symbol and algebraic operator for the 
simplest oscillatory process. As a symbolic form, i is an eigenform satisfying the equation 

i = 

One does not have an increment of time all alone as in classical t. One has it, a combination of 
an interval and the elemental dynamic that is time. With this understanding, we can return to the 
commutator for a discrete process and use it for the temporal increment. 

We found that discrete observation led to the commutator equation 

[x, Dx] = J(Ax) 2 /At 

which we will simplify to 

[q,p/m] = (Ax) 2 /At. 

taking q for the position x and p/m for velocity, the time derivative of position and ignoring the 
time shifting operator on the right hand side of the equation. 

Understanding that At should be replaced byiAt, and that, by comparison with the physics 
of a process at the Planck scale one can take 

(Ax) 2 /At = h/m, 

we have 

[q,p/m] = (Ax) 2 /iAt = —ih/m, 

whence 

[p, q\ = ih, 

and we have arrived at Heisenberg's fundamental relatiionship between position and momentum. 
This mode of arrival is predicated on the recognition that only it represents a true interval of 
time. In the notion of time there is an inherent clock or an inherent shift of phase that is making 
a synchrony in our ability to observe, a precise dynamic beneath the apparent dynamic of the 
observed process. Once this substitution is made, once the correct imaginary value is placed in 
the temporal circuit, the patterns of quantum mechanics appear. In this way, quantum mechanics 
can be seen to emerge from the discrete. 

The problem that we have examined in this section is the problem to understand the nature 
of quantum mechanics. In fact, we hope that the problem is seen to disappear the more we enter 
into the present viewpoint. A viewpoint is only on the periphery. The iterant from which the 
viewpoint emerges is in a superposition of indistinguishables, and can only be approached by 
varying the viewpoint until one is released from the particularities that a point of view contains. 
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3 Review of Non- Commutative Worlds 



Now we begin the introduction to non-commutative worlds and a general discrete calculus. Our 
approach begins in an algebraic framework that naturally contains the formalism of the calculus, 
but not its notions of limits or constructions of spaces with specific locations, points and trajecto- 
ries. Many patterns of physical law fit well into such an abstract framework. In this viewpoint one 
dispenses with continuum spacetime and replaces it by algebraic structure. Behind that structure, 
space stands ready to be constructed, by discrete derivatives and patterns of steps, or by start- 
ing with a discrete pattern in the form of a diagram, a network, a lattice, a knot, or a simplicial 
complex, and elaborating that structure until the specificity of spatio-temporal locations appear. 

Poisson brackets allow one to connect classical notions of location with the non-commutative 
algebra used herein. Below the level of the Poisson brackets is a treatment of processes and 
operators as though they were variables in the same context as the variables in the classical 
calculus. In different degrees one lets go of the notion of classical variables and yet retains 
their form, as one makes a descent into the discrete. The discrete world of non-commutative 
operators is a world linked to our familiar world of continuous and commutative variables. This 
linkage is traditionally exploited in quantum mechanics to make the transition from the classical 
to the quantum. One can make the journey in the other direction, from the discrete and non- 
commutative to the "classical" and commutative, but that journey requires powers of invention 
and ingenuity that are the subject of this exploration. It is our conviction that the world is basically 
simple. To find simplicity in the complex requires special attention and care. 

In starting from a discrete point of view one thinks of a sequence of states of the world 
S, S', S", S'", ■ ■ ■ where S' denotes the state succeeding S in discrete time. It is natural to suppose 
that there is some measure of difference DS^> = ,S( n+1 ) — S^ n \ and some way that states S and 
T might be combined to form a new state ST. We can thus think of world-states as operators 
in a non-commutative algebra with a temoporal derivative DS = S' — S. At this bare level 
of the formalism the derivative does not satisfy the Leibniz rule. In fact it is easy to verify 
that D(ST) = D(S)T + S'D{T). Remarkably, the Leibniz rule, and hence the formalisms of 
Newtonian calculus can be restored with the addition of one more operator J. In this instance J 
is a temporal shift operator with the property that S J = JS' for any state S. We then see that if 
VS = JD(S) = J(S' - S). then V(5T) = V(S)T + 5V(T) for any states S and T. In fact 
V(5) = JS' — JS = SJ — JS = [S, J], so that this adjusted derivative is a commutator in 
the general calculus of states. This, in a nutshell, is our approach to non-commutative worlds. 
We begin with a very general framework that is a non-numerical calculus of states and operators. 
It is then fascinating and a topic of research to see how physics and mathematics fit into the 
frameworks so constructed. 

A simplest and fundamental instance of these ideas is seen in the structure of i — We 
view i as an iterant (see Section 4), a discrete elementary dynamical system repeating in time 
the values {• • • — 1, +1, — 1, +1, • • • }. One can think of this system as resulting from the attempt 
to solve i 2 = —1 in the form i = Then one iterates the transformation x — > —1/x and 
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finds the oscillation from a starting value of +1 or —1. In this sense % is identical in concept to a 
primordial time. Furthermore the algebraic structure of the complex numbers emerges from two 
conjugate views of this discrete series as [—1, +1] and [+1,-1]. We introduce a temporal shift 
operator 7] such that rj[— 1, +1] = [+1, — l]r/ and rf — 1 (sufficient to this purpose). Then we can 
define i = [— 1, +l]r), endowing it with one view of the discrete oscillation and the sensitivity 
to shift the clock when interacting with itself or with another operator. See Sections 2 and 4 for 
the details of this reconstruction of the complex numbers. The point of the reconstruction for 
our purposes is that i becomes inextricably identified with elemental time, and so the physical 
substituion of it for t (Wick rotation) becomes, in this epistemology, an act of recognition of the 
nature of time. 

Constructions are performed in a Lie algebra A. One may take A to be a specific matrix Lie 
algebra, or abstract Lie algebra. If A is taken to be an abstract Lie algebra, then it is convenient 
to use the universal enveloping algebra so that the Lie product can be expressed as a commutator. 
In making general constructions of operators satisfying certain relations, it is understood that one 
can always begin with a free algebra and make a quotient algebra where the relations are satisfied. 

On A, a variant of calculus is built by defining derivations as commutators (or more generally 
as Lie products). For a fixed N in A one defines 

V N : A — > A 

by the formula 

\7 N F = [F, N] — FN — NF. 
Vat is a derivation satisfying the Leibniz rule. 

V N (FG) = V N (F)G + FV N (G). 

Discrete Derivatives are Replaced by Commutators. There are many motivations for replacing 
derivatives by commutators. If f(x) denotes (say) a function of a real variable x, and f(x) = 
f(x + h) for a fixed increment h, define the discrete derivative Df by the formula Df — (f — 
f) I h, and find that the Leibniz rule is not satisfied. One has the basic formula for the discrete 
derivative of a product: 

D(fg) = D(f)g + fD(g). 

Correct this deviation from the Leibniz rule by introducing a new non-commutative operator J 
with the property that 

fJ = Jf. 

Define a new discrete derivative in an extended non-commutative algebra by the formula 

V(/) = JD(f). 

It follows at once that 

V(/<7) = JD(f)g + JfD(g) = JD(f)g + fJD(g) = V(f)g + fV(g). 
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Note that 

V(/) = (Jf - Jf)/h = (fj - Jf)/h = [f, J/h}. 

In the extended algebra, discrete derivatives are represented by commutators, and satisfy the 
Leibniz rule. One can regard discrete calculus as a subset of non-commutative calculus based on 
commutators. 

Advanced Calculus and Hamiltonian Mechanics or Quantum Mechanics in a Non- Commutative 
World. In A there are as many derivations as there are elements of the algebra, and these deriva- 
tions behave quite wildly with respect to one another. If one takes the concept of curvature as the 
non-commutation of derivations, then A is a highly curved world indeed. Within A one can build 
a tame world of derivations that mimics the behaviour of flat coordinates in Euclidean space. 
The description of the structure of A with respect to these flat coordinates contains many of the 
equations and patterns of mathematical physics. 

The flat coordinates Q l satisfy the equations below with the Pj chosen to represent differentiation 
with respect to : 

[Q i ,Q i ] = o 
[Q\Pj\=S i5 . 

Here is the Kronecker delta, equal to 1 when i = j and equal to otherwise. Derivatives are 
represented by commutators. 

d l F = dF/dQ l = [F,Pl 

S, l F = dF/dP l = [Q\F]. 

Our choice of commutators guarantees that the derivative of a variable with respect to itself is 
one and that the derivative of a variable with respect to a distinct variable is zero. Furthermore, 
the commuting of the variables with one another guarantees that mixed partial derivatives are 
independent of the order of differentiation. This is a flat non-commutative world. 

Temporal derivative is represented by commutation with a special (Hamiltonian) element H 
of the algebra: 

dF/dt = [F, H] . 

(For quantum mechanics, take ihdA/dt = [A, H].) These non-commutative coordinates are the 
simplest flat set of coordinates for description of temporal phenomena in a non-commutative 
world. 

Hamilton's Equations are Part of the Mathematical Structure of Non-Commutative Ad- 
vanced Calculus. 

dPi/dt = [P h H] = -[H,Pi] = -dHjdQ 1 

dQ l /dt = [Q\H] = 8H/dP t . 

These are exactly Hamilton's equations of motion. The pattern of Hamilton's equations is built 
into the system. 



10 



The Simplest Time Series Leads to the Diffusion Constant and Heisenberg's Commuator. 

Consider a time series {Q, Q' , Q" , ...} with commuting scalar values. Let 

Q = VQ = JDQ = J(Q' - Q)/r 

where r is an elementary time step (If Q denotes a times series value at time t, then Q' denotes 
the value of the series at time t + r.). The shift operator J is defined by the equation Q J = JQ' 
where this refers to any point in the time series so that Q^ n >J = JQ^ n+1 ^ for any non-negative 
integer n. Moving J across a variable from left to right, corresponds to one tick of the clock. This 
discrete, non-commutative time derivative satisfies the Leibniz rule. 

This derivative V also fits a significant pattern of discrete observation. Consider the act of 
observing Q at a given time and the act of observing (or obtaining) DQ at a given time. Since Q 
and Q' are ingredients in computing (Q' — Q)/r, the numerical value associated with DQ, it is 
necessary to let the clock tick once, Thus, if one first observe Q and then obtains DQ, the result 
is different (for the Q measurement) if one first obtains DQ, and then observes Q. In the second 
case, one finds the value Q' instead of the value Q, due to the tick of the clock. 

1. Let QQ denote the sequence: observe Q, then obtain Q. 

2. Let QQ denote the sequence: obtain Q, then observe Q. 

The commutator [Q, Q] expresses the difference between these two orders of discrete mea- 
surement. In the simplest case, where the elements of the time series are commuting scalars, one 
has 

[Q, Q} = QQ-QQ = J(Q' - Qf /r. 
Thus one can interpret the equation 

[Q, Q] = Jk 

(k a constant scalar) as 

{Q'-Qf/r = k. 
This means that the process is a walk with spatial step 

A = ±Vk^ 

where k is a constant. In other words, one has the equation 

k = A 2 /r. 

This is the diffusion constant for a Brownian walk. A walk with spatial step size A and time step 
r will satisfy the commutator equation above exactly when the square of the spatial step divided 
by the time step remains constant. This shows that the diffusion constant of a Brownian process is 
a structural property of that process, independent of considerations of probability and continuum 
limits. 
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Thus we can write (ignoring the timeshift operator J) 

[Q,Q] = (AQ) 2 /r. 

If we work with physics at the Planck scale, then we can take r as the Planck time and AQ as the 
Planck length. Then 

(AQ) 2 /r = h/m 

where m is the Planck mass. However, we shall also Wick rotate the time from r to ir justifying 
ir on the principle (described above) that r should be multiplied by i to bring time into coinci- 
dence with an elemental time that is both a temporal operator (i) and a value (£). With this we 
obtain 

[Q, Q] = —ih/m 

or 

[mQ, Q] = ih, 

and taking P = mQ, we have finally 

[P, Q] = ih. 

Heisenberg's commutator for quantum mechanics is seen in the nexus of discrete physics and 
imaginary time. 

Schroedinger's Equation is Discrete. Here is how the Heisenberg form of Schroedinger's equa- 
tion fits in this context. Let J = (1 + ihHAt). Then = [tp, J / At], and we calculate 

= ip[(l + ihHAt) /At] - [(1 + ihH At) / At]ip = ih[i/j,H]. 

This is exactly the form of the Heisenberg equation. 

Dynamical Equations Generalize Gauge Theory and Curvature. One can take the general 
dynamical equation in the form 

dQ'/dt = Qi 

where {Q\, ■ ■ ■ , Qd} is a collection of elements of A. Write Qi relative to the flat coordinates via 
Qi = Pi — Ai. This is a definition of A± and dF/dQ' 1 = [F, PJ. The formalism of gauge theory 
appears naturally. In particular, if 

V l (F) = [F,G l ], 

then one has the curvature 

[V i ,V i ]F=[i2 ii ,F] 

and 

Rij = d{Aj — djAi + [Ai, Aj]. 

This is the well-known formula for the curvature of a gauge connection. Aspects of geometry 
arise naturally in this context, including the Levi-Civita connection (which is seen as a conse- 
quence of the Jacobi identity in an appropriate non-commutative world). 
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One can consider the consequences of the commutator [Q\ Qi] = gij, deriving that 

Qr = Q r + Frs Qs + T rst Q*Q\ 

where G r is the analogue of a scalar field, F rs is the analogue of a gauge field and T rst is the 
Levi-Civita connection associated with g^. This decompositon of the acceleration is uniquely 
determined by the given framework. 

Non-commutative Electromagnetism and Gauge Theory. One can use this context to revisit 
the Feynman-Dyson derivation of electromagnetism from commutator equations, showing that 
most of the derivation is independent of any choice of commutators, but highly dependent upon 
the choice of definitions of the derivatives involved. Without any assumptions about initial com- 
mutator equations, but taking the right (in some sense simplest) definitions of the derivatives 
one obtains a significant generalization of the result of Feynman-Dyson. We give this derivation 
in Section 5 of the present paper, using diagrammatic algebra to clarify the structure. In this 
derivation we use X to denote the position vector rather than Q, as above. 

Theorem With the appropriate [see below] definitions of the operators, and taking 

V 2 = df + d 2 2 +dl, H = X x X and E = d t X, one has 

1. X = E + X x H 

2. VH = 

3. d t H + V x E = H x H 

4. d t E-V xH = (<9 2 - V 2 )X 

The key to the proof of this Theorem is the definition of the time derivative. This definition 
is as follows 

d t F = F - EiXMF) = F — EiXi[F, X t ] 

for all elements or vectors of elements F. The definition creates a distinction between space and 
time in the non-commutative world. A calculation reveals that 

X = d t X + X x (X x X). 

This suggests taking E = d t X as the electric field, and B = X x X as the magnetic field so that 
the Lorentz force law 

X = E+X x B 

is satisfied. 

This result is applied to produce many discrete models of the Theorem. These models show that, 
just as the commutator [X, X] = Jk describes Brownian motion in one dimension, a generaliza- 
tion of electromagnetism describes the interaction of triples of time series in three dimensions. 
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Taking d t F = F — £jXj<9j (F) = F — EjXj [F, X»] as a definition of the partial derivative with 
respect to time is a natural move in this context because there is no time variable t in this non- 
commutative world. A formal move of this kind, matching a pattern from the commutative world 
to the mathematics of the non-commuative world is the theme of the Section 6 of this paper. In 
that section we consider the well known way to associate an operator to a product of commutative 
variables by taking a sum over all permutations of products of the operators corresponding to the 
individual variables. This provides a way to associate operator expressions with expressions in the 
commuative algebra, and hence to let a classical world correspond or map to a non-commutative 
world. To bind these worlds more closely, we can ask that the formulas for taking derivatives in 
the commuative world should have symmetrized operator product correspondences in the non- 
commuative world. In Section 6 we show how the resulting constraints are related to having a 
quadratic Hamiltonian (first order constraint) and to having a version of general relativity (31 
(second order constraint). Such constraints can be carried to all orders of derivatives, but the 
algebra of such constraints is, at the present time, in a very primitive state. We discuss some of 
the complexities of the constraint algebra in the Appendix to this paper. 

In Section 7 we discuss the relationship of the Bianchi identity in non-commutative worlds 
and its role in the classical derivation of Einstein's equations. This suggests other avenues of 
relationship between general relativity and non-commuative worlds. The reader may well aske at 
this point if we propose quantum gravity via this framework. In our judgement, it is too early to 
tell. 



Remark. While there is a large literature on non-commutative geometry, emanating from the 
idea of replacing a space by its ring of functions, work discussed herein is not written in that 
tradition. Non-commutative geometry does occur here, in the sense of geometry occuring in the 
context of non-commutative algebra. Derivations are represented by commutators. There are 
relationships between the present work and the traditional non-commutative geometry, but that 
is a subject for further exploration. In no way is this paper intended to be an introduction to 
that subject. The present summary is based on lfT5l [TTl [T%1 fT9l l20l |2T1 |22l l23l [24l [251l and the 
references cited therein. 



The following references in relation to non-commutative calculus are useful in comparing 
with the present approach [|2l[51|71|29]|. Much of the present work is the fruit of a long series of 
discussions with Pierre Noyes. paper [|27l also works with minimal coupling for the Feynman- 
Dyson derivation. The first remark about the minimal coupling occurs in the original paper by 
Dyson fl], in the context of Poisson brackets. The paper [8] is worth reading as a companion to 
Dyson. It is the purpose of this summary to indicate how non-commutative calculus can be used 
in foundations. 
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4 Iterants, Discrete Processes and Matrix Algebra 



The primitive idea behind an iterant is a periodic time series or "waveform" 

• • • abababababab ■ ■ ■ . 

The elements of the waveform can be any mathematically or empirically well-defined objects. 
We can regard the ordered pairs [a, b] and [b, a] as abbreviations for the waveform or as two 
points of view about the waveform (a first or b first). Call [a, b] an iterant. One has the collection 
of transformations of the form T[a, b] = [ka, k _1 b] leaving the product ab invariant. This tiny 
model contains the seeds of special relativity, and the iterants contain the seeds of general matrix 
algebra! Since this paper has been a combination of discussions of non-commutativity and time 
series, we include this appendix on iterants. A more complete discussion will appear elsewhere. 
For related discussion see |E1 HSl HB HH US I2E1 - 

Define products and sums of iterants as follows 

[a, b] [c, d] = [ac, bd] 

and 

[a, b] + [c, d] = [a + c,b + d]. 

The operation of juxtapostion is multiplication while + denotes ordinary addition in a category 
appropriate to these entities. These operations are natural with respect to the structural juxtapo- 
sition of iterants: 

...abababababab... 
...cdcdcdcdcdcd... 

Structures combine at the points where they correspond. Waveforms combine at the times where 
they correspond. Iterants conmbine in juxtaposition. 

If • denotes any form of binary compositon for the ingredients (a,b,...) of iterants, then we 
can extend • to the iterants themselves by the definition [a, b] • [c, d] = [a»c,b»d]. In this section 
we shall first apply this idea to Lorentz transformations, and then generalize it to other contexts. 

So, to work: We have 

[t-x,t + x] = [t,t] + [-x,x] = t[l, 1] 
Since [1, l][a, b] = [la, lb] = [a, b] and [0, 0}[a, b] = [0, 0], we shall write 

1 = [1,1] 

and 

0= [0,0]. 
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Let 

o =[-1,1]- 

a is a significant iterant that we shall refer to as a polarity. Note that 

aa = 1. 

Note also that 

[t — X, t + x] = t + £0". 

Thus the points of spacetime form an algebra analogous to the complex numbers whose elements 
are of the form t + xa with aa = 1 so that 

(t + xa)(t' + x'a) = tif + xx' + (tx 1 + xt')a. 

In the case of the Lorentz transformation it is easy to see the elements of the form [k, k" 1 } translate 
into elements of the form 

T(v) = [(1 + v)/yfty^P), (1 - v)/y^T^)} = [k, k- 1 }. 

Further analysis shows that v is the relative velocity of the two reference frames in the physical 
context. Multiplication now yields the usual form of the Lorentz transform 

T k (t + xa) = T{v)(t + xa) 

= (1/7(1 -u 2 ) - vcx/y/(l-v*))(t + xa) 
= (t- xv)/^{l-v 2 ) + (x - vt)a/^{l-v 2 ) 
= t + x a. 

The algebra that underlies this iterant presentation of special relativity is a relative of the 
complex numbers with a special element a of square one rather than minus one (i 2 = — 1 in the 
complex numbers). 

The appearance of a square root of minus one unfolds naturally from iterant considerations. 
Define the "shift" operator 77 on iterants by the equation 

rj[a,b] = [b,a]rj 

with rj 2 = 1. Sometimes it is convenient to think of 7] as a delay opeator, since it shifts the 
waveform ...ababab... by one internal time step. Now define 

i=[-l,l]v 

We see at once that 

a = [-1, 1M-1, % = [-1, i][i, -% 2 = [-1, i][i, -1] = [-1, -1] = -1. 
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Thus 



ii = — 1. 

Here we have described i in a new way as the superposition of the waveform e = [—1,1] and the 
temporal shift operator rj. By writing i = er/ we recognize an active version of the waveform that 
shifts temporally when it is observed. This theme of including the result of time in observations 
of a discrete system occurs at the foundation of our construction. 



4.1 MATRIX ALGEBRA VIA ITERANTS 

Matrix algebra has some strange wisdom built into its very bones. Consider a two dimensional 
periodic pattern or "waveform." 

. . . abababababababab. . . 
. . . cdcdcdcdcdcdcdcd. . . 
. . . abababababababab. . . 
. . . cdcdcdcdcdcdcdcd. . . 
. . . abababababababab. . . 



a b \ f b a\ f c d\ ( d c 
c d)\d c )\a b )\b a 

Above are some of the matrices apparent in this array. Compare the matrix with the "two dimen- 
sional waveform" shown above. A given matrix freezes out a way to view the infinite waveform. 
In order to keep track of this patterning, lets write 



[a, d] + [b, c}r] 

where 

fay] 

and 

V = 



a b 
c d 



x 
y 

1 

1 



The four matrices that can be framed in the two-dimensional wave form are all obtained from the 
two iterants [a, d] and [b, c] via the shift operation r)[x,y] = [y, x]r) which we shall denote by an 
overbar as shown below 

[x,y] = [y,x]. 
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Letting A = [a, d] and B = [b, c], we see that the four matrices seen in the grid are 

A + Brj, B + At), B + Ar), A + Br). 

The operator 77 has the effect of rotating an iterant by ninety degrees in the formal plane. Ordinary 
matrix multiplication can be written in a concise form using the following rules: 

1)7] = 1 

T)Q = Qt] 

where Q is any two element iterant. 

For example, let e = [—1, 1] so that e = — e and ee = [1, 1] = 1. Let 

i = erj. 

Then 

ii = er]er] = eer/r] = e(— e) = — ee = —1. 
We have reconstructed the square root of minus one in the form of the matrix 

i = er] = [—1, l]r] = 

More generally, we see that 

(A + Brj)(C + Dr)) = (AC + BD) + (AD + BC)r) 
writing the 2x2 matrix algebra as a system of hypercomplex numbers. Note that 

(A + Brj)(A - Brj) = A~A — B~B 

The formula on the right corresponds to the determinant of the matrix. Thus we define the 
conjugate of A + Br/ by the formula 

A + Br] = A- Brj. 
These patterns generalize to higher dimensional matrix algebra. 

It is worth pointing out the first precursor to the quaternions: This precursor is the system 

{±l,±e, ±7j, ±i}. 

Here ee = 1 = 7777 while i = er] so that ii = —1. The basic operations in this algebra are those of 
epsilon and eta. Eta is the delay shift operator that reverses the components of the iterant. Epsilon 
negates one of the components, and leaves the order unchanged. The quaternions arise directly 
from these two operations once we construct an extra square root of minus one that commutes 
with them. Call this extra root of minus one Then the quaternions are generated by 

{i = eT j J = yf-[e, k = V^lr)} 
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with 



i 2 = j 2 = k 2 = ijk = — 1. 



The "right" way to generate the quaternions is to start at the bottom iterant level with boolean 
values of and 1 and the operation EXOR (exclusive or). Build iterants on this, and matrix 
algebra from these iterants. This gives the square root of negation. Now take pairs of values from 
this new algebra and build 2x2 matrices again. The coefficients include square roots of negation 
that commute with constructions at the next level and so quaternions appear in the third level of 
this hierarchy. 



4.2 Matrix Algebra in General 

Construction of matrix algebra in general proceeds as follows. Let M be an n x n matrix over a 
ring R. Let M = (m^) denote the matrix entries. Let ir be an element of the symmetric group 
S n so that 7Ti, 7T2, ■ ■ • , 7r n is a permuation of 1, 2, ■ • • , n. Let v — i>2, • • ■ , v n ) denote a vector 
with these components. Let A(v) denote the diagonal matrix whose i — th diagonal entry is Uj. 
Let v w = (vjtj, ■ ■ ■ , v nn ). Let A n (v) = A(v n ). Let A denote any diagonal matrix and A n denote 
the corresponding permuted diagonal matrix as just described. Let [it] denote the permutation 
matrix obtained by taking the i — th row of [ir] to be the 7Tj — th row of the identity matrix. Note 
that [it] A = A n [it] . For each element it of S n define the vector v (M, it) = (m kl , • • • , m n7Tn ) and 
the diagonal matrix A[M] W = A(v(M, it)). 



Theorem. M = (l/(n - l)!)E 7r65n A[M] 7r [7r]. 



The proof of this theorem is omitted here. Note that the theorem expresses any square matrix 
as a sum of products of diagonal matrices and permutation matrices. Diagonal matrices add 
and multiply by adding and multiplying their corresponding entries. They are acted upon by 
permutations as described above. This means that any matrix algebra can be embedded in an 
algebra that has the structure of a group ring of the permutation group with coefficients A in an 
algebra (here the diagonal matrices) that are acted upon by the permutation group, and following 
the rule [tt]A = A 77 [tt] . This is a full generalization of the case n = 2 described in the last section. 



For example, we have the following expansion of a 3 x 3 matrix: 



a 


b 






( a 








d 


e 


f 





e 





9 


h 


I) 


Vo 





k 










c 


d 











h 












c 





e 





9 









6 
d 

A; 



a 
0/ 
h 
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Here, each term factors as a diagonal matrix multiplied by a permutation matrix as in 

a \ /a0 0\/l0 0\ 

oo/ = o/ o ooi . 

h / \0 0/i/\0 10/ 

It is amusing to note that this theorem tells us that up to the factor of l/(n — 1)! a unitary matrix 
that has unit complex numbers as its entries is a sum of simpler unitary transformations factored 
into diagonal and permutation matrices. In quantum computing parlance, such a unitary matrix is 
a sum of products of phase gates and products of swap gates (since each permutation is a product 
of transpositions). 

A reason for discussing these formulations of matrix algebra in the present context is that one 
sees that matrix algebra is generated by the simple operations of juxtaposed addition and multi- 
plication, and by the use of permutations as operators. These are unavoidable discrete elements, 
and so the operations of matrix algebra can be motivated on the basis of discrete physical ideas 
and non-commutativity. The richness of continuum formulations, infinite matrix algebra, and 
symmetry grows naturally out of finite matrix algebra and hence out of the discrete. 



5 Generalized Feynman Dyson Derivation 

In this section we assume that specific time- varying coordinate elements Xi,X 2 , X 3 of the alge- 
bra A are given. We do not assume any commutation relations about X±, X 2 , X%. 

In this section we no longer avail ourselves of the commutation relations that are in back of 
the original Feynman-Dyson derivation. We do take the definitions of the derivations from that 
previous context. Surprisingly, the result is very similar to the one of Feynman and Dyson, as we 
shall see. 

Here A x Bis the non-commutative vector cross product: 

(A x B) k = El^jkAiBj. 

(We will drop this summation sign for vector cross products from now on.) Then, with B = 
Ixl,we have 

Bk = tijkXiXj = (l/2)eijk[Xi, Xj\. 

The epsilon tensor is defined for the indices {i, j, k} ranging from 1 to 3, and is equal to if 
there is a repeated index and is ortherwise equal to the sign of the permutation of 123 given by 
ijk. We represent dot products and cross products in diagrammatic tensor notation as indicated 
in Figure 2 and Figure 3. In Figure 2 we indicate the epsilon tensor by a trivalent vertex. The 
indices of the tensor correspond to labels for the three edges that impinge on the vertex. The 
diagram is drawn in the plane, and is well-defined since the epsilon tensor is invariant under 
cyclic permutation of its indices. 
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We will define the fields E and B by the equations 

B = X x X and E = d t X. 

We will see that E and B obey a generalization of the Maxwell Equations, and that this gen- 
eralization describes specific discrete models. The reader should note that this means that a 
significant part of the form of electromagnetism is the consequence of choosing three coordinates 
of space, and the definitions of spatial and temporal derivatives with respect to them. The back- 
ground process that is being described is otherwise aribitrary, and yet appears to obey physical 
laws once these choices are made. 

In this section we will use diagrammatic matrix methods to carry out the mathematics. In 
general, in a diagram for matrix or tensor composition, we sum over all indices labeling any edge 
in the diagram that has no free ends. Thus matrix multiplication corresponds to the connecting of 
edges between diagrams, and to the summation over common indices. With this interpretation of 
compositions, view the first identity in Figure 2. This is a fundmental identity about the epsilon, 
and corresponds to the following lemma. 




ABC ABC ABC 




Figure 2: Epsilon Identity 

Lemma. (View Figure 2) Let e ijfc be the epsilon tensor taking values 0, 1 and —1 as follows: 
When ijk is a permuation of 123, then e^. is equal to the sign of the permutation. When ijk 
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contains a repetition from {1,2, 3}, then the value of epsilon is zero. Then e satisfies the labeled 
identity in Figure 2 in terms of the Kronecker delta. 

The proof of this identity is left to the reader. The identity itself will be referred to as the epsilon 
identity. The epsilon identity is a key structure in the work of this section, and indeed in all 
formulas involving the vector cross product. 

The reader should compare the formula in this Lemma with the diagrams in Figure 2. The 
first two diagram are two versions of the Lemma. In the third diagram the labels are capitalized 
and refer to vectors A, B and C. We then see that the epsilon identity becomes the formula 

A x (B x C) = {A • C)B - {A • B)C 

for vectors in three-dimensional space (with commuting coordinates, and a generalization of this 
identity to our non-commutative context. Refer to Figure 3 for the diagrammatic definitions of 
dot and cross product of vectors. We take these definitions (with implicit order of multiplication) 
in the non-commutative context. 



A.B = A B 


-a.F =[f, 






V 


Ax B = A B 


Y 


f =<a t F + 

/ / 


X[F,X] 











VxF = ^ F 




= [F,X] =-[F,X] 




Figure 3: Defining Derivatives 



Remarks on the Derivatives. 
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1. Since we do not assume that [X;,X,] = nor do we assume [Xi,Xj] = 0, it will not 
follow that E and B commute with the X { . 

2. We define 

di(F) = [F,Xi], 

and the reader should note that, these spatial derivations are no longer flat in the sense of 
section 1 (nor were they in the original Feynman-Dyson derivation). See Figure 3 for the 
diagrammatic version of this definition. 

3. We define d t = d/dtby the equation 

d t F = F- EiXidi(F) — F — £ 4 X 4 [F, X t ] 

for all elements or vectors of elements F. We take this equation as the global definition of 
the temporal partial derivative, even for elements that are not commuting with the Xj. This 
notion of temporal partial derivative d t is a least relation that we can write to describe the 
temporal relationship of an arbitrary non-commutative vector F and the non-commutative 
coordinate vector X. See Figure 3 for the diagrammatic version of this definition. 

4. In defining 

d t F = F-Z i X i d i (F), 

we are using the definition itself to obtain a notion of the variation of F with respect to time. 
The definition itself creates a distinction between space and time in the non-commutative 
world. 

5. The reader will have no difficulty verifying the following formula: 

d t (FG) = d t (F)G + Fd t (G) + ZAiFMG). 

This formula shows that d t does not satisfy the Leibniz rule in our non-commutative con- 
text. This is true for the original Feynman-Dyson context, and for our generalization of 
it. All derivations in this theory that are defined directly as commutators do satisfy the 
Leibniz rule. Thus d t is an operator in our theory that does not have a representation as a 
commutator. 

6. We define divergence and curl by the equations 

VmB = Y* =l d i {B i ) 

and 

(V x E) k = e ljk d i (E j ). 
See Figure 3 and Figure 5 for the diagrammatic versions of curl and divergence. 
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Now view Figure 4. We see from this Figure that it follows directly from the definition of the 
time derivatives (as discussed above) that 

X = d t X + X x(X x X). 

This is our motivation for defining 

E = d t X 

and 

B = X xX. 

With these definition in place we have 

X = E + X x B, 

giving an analog of the Lorentz force law for this theory. 

Just for the record, look at the following algebraic calculation for this derivative: 

F = a t F + E i X i [F,X i ] 

= d t F + V i (X i FX i -X i X i F) 
= d t F + ^{XiFXi - XiFX) + X.FiX - XiX.F 

Hence 

F = d t F + X x F + (X • F)X - (X • X)F 
(using the epsilon identity). Thus we have 

X = d t X + X x (X x X) + (X • X)X - (X • X)X, 

whence 

X = d t X + X x(X x X). 
In Figure 5, we give the derivation that B has zero divergence. 

Figures 6 and 7 compute derivatives of B and the Curl of E, culminating in the formula 

d t B + V x E = B x B. 

In classical electromagnetism, there is no term B x B. This term is an artifact of our non- 
commutative context. In discrete models, as we shall see at the end of this section, there is 
no escaping the effects of this term. 

Finally, Figure 8 gives the diagrammatic proof that 

d t E-V x B = (d 2 t - V 2 )X. 

This completes the proof of the Theorem below. 

24 



F = ^ t F + X[F,X] 

/ / V 

F = c) t F + X FX - XX F 

x = / r) t x + xxx - xxx 

• • • • 

= Q.X + X X X 

; y 

X* = 'StX + X x (X x X) 
Figure 4: The Formula for Acceleration 

Electromagnetic Theorem With the above definitions of the operators, and taking 

V 2 = dj + d\ + <9 2 , B = lxi and E = d t X we have 

1. X = E + X x B 

2. VB = 

3. d t B + VxE = BxB 

4. d t E-V x B = (9 2 - V 2 )X 

Remark. Note that this Theorem is a non-trivial generalization of the Feynman-Dyson derivation 
of electromagnetic equations. In the Feynman-Dyson case, one assumes that the commutation 
relations 

[X i ,X j ]=Q 

and 

[X h Xj] = 5ij 

are given, and that the principle of commutativity is assumed, so that if A and B commute with 
the Xi then A and B commute with each other. One then can interpret di as a standard derivative 
with di(Xj) = Sij. Furthermore, one can verify that Ej and Bj both commute with the Xi. From 
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E = ^ t X B = X x X 

X = E + Xx B 







VB = [B , X] 




LJ 




• • • • • 


• • • 


=BX - XB=XXX 


- XXX = 






V»B = 





Figure 5: Divergence of B 



this it follows that d t (E) and d t (B) have standard intepretations and that B x B = 0. The above 
formulation of the Theorem adds the description of E as d t (X), a non-standard use of d t in the 
original context of Feyman-Dyson, where d t would only be defined for those A that commute 
with Xi. In the same vein, the last formula d t E - Vx5 = (<9 t 2 — V 2 )X gives a way to express 
the remaining Maxwell Equation in the Feynman-Dyson context. 

Remark. Note the role played by the epsilon tensor throughout the construction of gener- 
alized electromagnetism in this section. The epsilon tensor is the structure constant for the Lie 
algebra of the rotation group SO (3). If we replace the epsilon tensor by a structure constant 
fijk for a Lie algebra Qof dimension d such that the tensor is invariant under cyclic permutation 
(fijk = fkijX then most of the work in this section will go over to that context. We would then 
have d operator/variables X 1: - ■ ■ X d and a generalized cross product defined on vectors of length 
d by the equation 

(AxB) k = fijkAiBj. 
The Jacobi identity for the Lie algebra Q implies that this cross product will satisfy 

A x [B x C) = (A x B) x C + [B x (A] x C) 

where 

{[B X {A} X C) r = fklrfijkABkCj. 

This extension of the Jacobi identity holds as well for the case of non-commutative cross product 
defined by the epsilon tensor. It is therefore of interest to explore the structure of generalized 
non-commutative electromagnetism over other Lie algebras (in the above sense). This will be the 
subject of another paper. 
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'(DtB = B + X [X, B] 



B = (1/2)[X,X* = [*X , X ] 



= [E,X] + [ XxB, X ] 



= - yx E + [ X B , X ] 




Figure 6: Computing B 



5.1 Discrete Thoughts 

In the hypotheses of the Electromagnetic Theorem, we are free to take any non-commutative 
world, and the Electromagnetic Theorem will satisfied in that world. For example, we can take 
each Xi to be an arbitary time series of real or complex numbers, or bitstrings of zeroes and ones. 
The global time derivative is defined by 

F = J(F f -F) = [F,J], 

where F J = JF'. This is the non-commutative discrete context discussed in sections 2 and 3. 
We will write 

F = JA(F) 

where A(F) denotes the classical discrete derivative 

A(F) — F' — F. 

With this interpretation X is a vector with three real or complex coordinates at each time, and 

B = X xX = J 2 A(X') x A(X) 

while 

E = X-Xx{XxX) = J 2 A 2 {X) - J 3 A(X") x (A(X') x A(X)). 

Note how the non-commutative vector cross products are composed through time shifts in this 
context of temporal sequences of scalars. The advantage of the generalization now becomes 
apparent. We can create very simple models of generalized electromagnetism with only the 
simplest of discrete materials. In the case of the model in terms of triples of time series, the 
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^B +V* E=X[X,B] + [XB,X] 




= -XXB + XXB ( Note that X B = B X ) 



= X X B =BxB 




'E) t B +Vx E = Bx B 



Figure 7: Curl of E 



generalized electromagnetic theory is a theory of measurements of the time series whose key 
quantities are 

A(X') x A{X) 

and 

A(X") x (A(X') x A(X)). 

It is worth noting the forms of the basic derivations in this model. We have, assuming that F 
is a commuting scalar (or vector of scalars) and taking Aj = X[ — X i7 

d t (F) = [F t Xi] = [F, J A,} = FJAi - JA,F = J{F'A i - A,F) = FA, 

and for the temporal derivative we have 

d t F = J[l - J A' • A]A(F) 

where A = (A l5 A 2 , A 3 ). 



6 Constraints - Classical Physics and General Relativity 

The program here is to investigate restrictions in a non-commutative world that are imposed by 
asking for a specific correspondence between classical variables acting in the usual context of 
continuum calculus, and non-commutative operators corresponding to these classical variables. 
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E = t X ^^E^fx 

VxB= ft XX 

Y 

XX t'c) XX 

2 2 # 

^E - VxB = CQ t - v )x 



Figure 8: Curl of £ 



By asking for the simplest constraints we find the need for a quadratic Hamiltonian and a remark- 
able relationship with Einstein's equations for general relativity 01 HJ. There is a hierarchy of 
constraints of which we only analyze the first two levels. An appendix to this paper indicates a 
direction for exploring the algebra of the higher constraints. 

If, for example, we let x and y be classical variables and X and Y the corresponding non- 
commutative operators, then we ask that x n correspond to X n and that y n correspond to Y n for 
positive integers n. We further ask that linear combinations of classical variables correspond to 
linear combinations of the corresponding operators. These restrictions tell us what happens to 
products. For example, we have classically that (x + y) 2 = x 2 + 2xy + y 2 . This, in turn must 
correspond to (X + Y) 2 = X 2 + XY + YX + Y 2 . From this it follows that 2xy corresponds to 
XY + YX. Hence xy corresponds to 

{XY} = (XY + YX)/2. 

By a similar calculation, if x±, X2, ■ ■ • ,x n are classical variables, then the product 2ix 2 ■ • -x n 
corresponds to 

{X 1 X 2 ■■■X n } = {l/n^ eSn X ai X a2 ■■■X an . 

where S n denotes all permutations of 1, 2, ■ ■ • ,n. Note that we use curly brackets for these sym- 
metrizers and square brackets for commutators as in [A, B] = AB — BA. 

We can formulate constraints in the non-commutative world by asking for a correspondence 
between familiar differentiation formulas in continuum calculus and the corresponding formulas 
in the non-commutative calculus, where all derivatives are expressed via commutators. We will 
detail how this constraint algebra works in the first few cases. Exploration of these constraints 
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has been pioneered by Tony Deakin [310. The author of this paper and Tony Deakin are planning 
a comprehensive paper on the consequences of these contraints in the interface between classical 
and quantum mechanics. 

Recall that the temporal derivative in a non-commutative world is represented by commutator 
with an operator H that can be intrepreted as the Hamiltonian operator in certain contexts. 

Q=[Q,H]. 

For this discussion, we shall take a collection Q 1 , Q 2 , ■ ■ ■ , Q n of operators to represent spatial 
coordinates q 1 , q 2 , ■ ■ ■ , q n . The Q % commute with one another, and the derivatives with respect to 
Q 1 are represented by operators P % so that 

We also write 

de/dP i = e* = [Q*,e]. 

To this purpose, we assume that [Q l , P J ] = and that the P J commute with one another (so 
that mixed partial derivatives with respect to the Q l are independent of order of differentiation). 

Note that 

Qi = [Q\ H) = H\ 

It will be convenient for us to write W in place of Q l in the calculations to follow. 
The First Constraint. The first constraint is the equation 

e = = {H i e t }. 

This equation expresses the symmetrized version of the usual calculus formula 9 = q i 9 i . It is 
worth noting that the first constraint is satisfied by the quadratic Hamiltonian 

H=±{g ij P i Pi + P i P3g ij ) 

where g t j = gji and the commute with the Q k . We leave the verification of this point to the 
reader, and note that the fact that the quadratic Hamiltonian does satisfy the first constraint shows 
how the constraints bind properties of classical physics (in this case Hamiltonian mechanics) to 
the non-commutative world. 

The Second Constraint. The second constraint is the symmetrized analog of the second tempo- 
ral derivative: 

e = {we t } + {u'hh-),,}. 

However, by differentiating the first constraint we have 

e = {WQt} + {//'{//•'(-)„}} 

Thus the second constraint is equivalent to the equation 

{H'iWQij}} = {//'//•' (-), 7 ), 
We now reformulate this version of the constraint in the following theorem. 
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Theorem. The second constraint in the form {H l {H j Qij}} = {H l H j Qij} is equivalent to the 
equation 

[[Q ij ,Hl],H i ] = 0. 

Proof. We can shortcut the calculations involved in proving this Theorem by looking at the 
properties of symbols A,B,C such that AB = BA, ACB = BCA. Formally these mimic 
the behaviour of A = H\ B = W , C = in the expressions H l WQij and WQijffl since 
Qij = Qji, and the Einstein summation convention is in place. Then 

{A{BC}} = ^{A{BC + CB) + (BC + CB)A) 

= ^(ABC + ACB + BCA + CBA), 

{ABC} = -(ABC + ACB + BAC + BCA + CAB + CBA). 
6 

So 

{ABC} - {A{BC}} = ^(-ABC - ACB + 2BAC - BCA + 2CAB - CBA) 

= -^(ABC - 2 ACB + CAB) 
= jf;(ABC - 2BCA + CBA) 
= ^{A{BC - CB) + (CB - BC)A) 
= ±(A[B,C]-[B,C]A) 
= ^[A,[B,C]}. 

Thus the second constraint is equivalent to the equation 

[H i ,[H\e ij ]]=0. 
This in turn is equivalent to the equation 

[[Q ij ,H*],H i ]=0, 

completing the proof of the Theorem. 
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Remark. If we define 

v l (e) = [e,^] = [e,g i ] 

then this is the natural covariant derivative that was described in the introduction to this paper. 
Thus the second order constraint is 

V'(V ((-),,) = 0. 

Note that 

= v i (v i (e ii ) + [[fr < , J ffJ'],e y ] 
= v i (v i (e ii ) + [[^ i ,^'],e y ]. 

Hence the second order constraint is equivalent to the equation 

[[H i ,H*),e ij \ = 0. 
This equation weaves the curvature of V with the flat derivatives of 0. 
A Relationship with General Relativity. Again, if we define 

w i (e) = [e,H i ] = [e,Qi] 

then this is the natural covariant derivative that was described in the introduction to this paper. 
Thus the second order constraint is 

v*(v j '(e y ) = o. 

If we use the quadratic Hamiltonian H = j(gijP l P j + P' l P j g i j) as above, then with = g lm 
the second constraint becomes the equation 

9 (9 9,jku)v - u. 

Deakin and Kilmister [0]| observe that this last equation specializes to a fourth order version of 
Einstein's field equation for vacuum general relativity: 

, 2 

K ab = g ef (R ab ;ef + -RaeRfb) = 

where a,b,e,f = 1, 2, ■ ■ • n and R is the curvature tensor corresponding to the metric g a f,. This 
equation has been studied by Deakin in [4J. It remains to be seen what the consequences for 
general relativity are in relation to this formulation, and it remains to be seen what the further 
consequences of higher order constraints will be. 

The algebra of the higher order constraints is under investigation at this time. 



32 



7 On the Algebra of Constraints 



We have the usual advanced calculus formula 9 = q i 6 i . We shall define h? = q i so that we can 
write 9 = h % 6i. We can then calculate successive derivatives with denoting the n-th temporal 
derivative of 9. 

9 W = h l 9i 
0( 2 ) = h^di + WhPdij 

0(3) = h i(2) e , + 3h m h 3Q.. + h i h j h% jk 

The equality of mixed partial derivatives in these calculations makes it evident that one can use a 
formalism that hides all the superscripts and subscripts (i, j, k, ■ ■ ■ )). In that simplified formalism, 
we can write 

9 {1) = h9 
0( 2 ) = h^9 + h 2 9 

0(3) = h (2) Q + 3h (l) h Q + tfg 

e& = h A 9 + Qh 2 9h^ + 39h^ 2 + Ah9h^ + ^( 3 ) 

Each successive row is obained from the previous row by applying the identity = h9 in 
conjunction with the product rule for the derivative. 

This procedure can be automated so that one can obtain the formulas for higher order deriva- 
tives as far as one desires. These can then be converted into the non-commutative constraint 
algebra and the consequences examined. Further analysis of this kind will be done in a sequel to 
this paper. 

The interested reader may enjoy seeing how this formalism can be carried out. Below we il- 
lustrate a calculation using Mathematical M , where the program already knows how to formally 
differentiate using the product rule and so only needs to be told that 6^> = h9. This is said in 
the equation T'[x] = H[x]T[x] where T[x] stands of 9 and H[x] stands for h with x a dummy 
variable for the differentiation. Here D[T[a;], x] denotes the derivative of T[x] with respect to x, 
as does T'[x], 

In the calculation below we have indicated five levels of derivative. The structure of the 
coefficients in this recursion is interesting and complex territory. For example, the coefficients of 
-ff[x] n T[x]-fT'[x] = h n 9h' are the triangular numbers {1, 3, 6, 10, 15, 21, • • • } but the next series 
are the coefficients of H[x] n T[x]H' [x] 2 = h n 9h' 2 , and these form the series 

{1, 3, 15, 45, 105, 210, 378, 630, 990, 1485, 2145, • • • }. 

This series is eventually constant after four discrete differentiations. This is the next simplest 
series that occurs in this structure after the triangular numbers. To penetrate the full algebra 
of constraints we need to understand the structure of these derivatives and their corresponding 
non-commutative symmetrizations. 
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T'[x]:=H[x]T[x] 

D[T[x],x] 
D[D[T[x] 

D[D[D[T[x],x],x],x] 
£)[D[D[Z)[r[x], x], x], x], x] 
D[D[D[D[D[T[x] 

ff[x]T[x] 

#[x] 2 T[x] + T[x]#'[x] 

Jf[x] 3 T[x] + 3#[x]T[x]if'[x] + T[x]H"[x] 

H[x] 4 T[x] + 6fl'[a;] 2 r[a;] J ff / [ar] + 3T[x]H'[x] 2 + 4i7[a;]T[a;]/f"[a;] + T[x]H^[x] 

H [x] 5 T [x] + 10H [x] 3 T [x] H' [x] + 15H[x]T [x] H' [x] 2 + 10 H[x} 2 T [x] H" [x] + 1 OT [x] H' [x] H" [x] + 
5H[x]T[x]H®[x] + T[x}H^[x] 

7.1 Algebra of Constraints 

In this section we work with the hidden index conventions described before in the paper. In this 
form, the classical versions of the first two constraint equations are 

1. 9 = 9h 

2. = 6h 2 + 6h 

In order to obtain the non-commutative versions of these equations, we replace h by H and 
9 by where the capitalized versions are non-commuting operators. The first and second con- 
straints then become 

1. {0} = {QH} = ~(QH + HQ) 

2. {0} = {QH 2 } + {QH} = \(QH 2 + HQH + H 2 Q) + \{QH + if 6) 
Proposition. The Second Constraint is equivalent to the commutator equation 

[[Q,H},H] = 0. 

Proof. We identify 

{©}' = {0} 

and 

{Q}' = {{QH}H} + {QH}. 



34 



So we need 

{QH 2 } = {{QH}H}. 
The explicit formula for {{QH} H} is 

{{QH}H} = ^({QH}H + H{QH}) = ^(9HH + HQH + HQH + HHQ). 

Thus we require that 

\{QH 2 + HQH + H 2 Q) = ^(9HH + HQH + HQH + HHQ). 

which is equivalent to 

QH 2 + H 2 Q - 2HQH = 0. 

We then note that 

[[0, H},H} = (QH - HQ)H - H(QH - HQ) = QH 2 + H 2 Q - 2HQH. 
Thus the final form of the second constraint is the equation 

[[Q,H],H]=0.// 

The Third Constraint. We now go on to an analysis of the third constraint. The third constraint 
consists in the the two equations 

1. {0} = {QH 3 } + 3{QHH} + {QH} 

2. {0} = {0}« where 

{0}« = {{QH}H 2 } + 2{QHH} + {{QH}H} + 

Proposition. The Third Constraint is equivalent to the commutator equation 

[H 2 ,[H,Q]] = [H, [H,Q]]-2[H, [H, ©]]. 
Proof. We demand that {©} = {©}* and this becomes the longer equation 

{QH 3 } + 3{QHH} + {QH} = {{QH}H 2 } + 2{QHH} + {{QH}H} + {QH} 
This is equivalent to the equation 

{QH 3 } + {QHH} = {{QH}H 2 } + {{QH}H} 
This, in turn is equivalent to 

{QH 3 } - {{QH}H 2 } = {{QH}H} - {QHH} 
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This is equivalent to 

(1/A)(H 3 Q+H 2 QH+HQH 2 +QH 3 )-(1/Q)(H 2 (HQ+QH)+H(HQ+QH)H+(HQ+QH)H 2 ) 

= (1/2)(H(1/2)(HQ+QH) + (1/2)(HQ+QH)H)-(1/6)(HHQ+HQH+HHQ+HQH+QHH+QHH) 
This is equivalent to 

3(H 3 Q + H 2 OH + HQH 2 + QH 3 ) - 2(H 3 Q + 2H 2 QH + 2HQH 2 + QH 3 ) 

= 3(HHQ + HQH + HQH + QHH) - 2(HHQ + HQH + HHQ + HQH + QHH + QHH) 
This is equivalent to 

H 3 Q - H 2 QH - HQH 2 + QH 3 
= (HHQ + HQH + HQH + QHH) - 2(HHQ + QHH) 
The reader can now easily verify that 

[H 2 , [H, Q]] = H 3 Q - H 2 QH - HQH 2 + QH 3 

and that 

[H, [H, 0]] - 2[H, [H, Q}} = (HHQ + HQH + HQH + QHH) - 2(HHQ + QHH) 
Thus we have proved that the third constraint equations are equivalent to the commutator equation 

[H 2 ,[H,Q]} = [H, [H,Q]]-2[H, [H, 0]] 
This completes the proof of the Proposition. // 

Discussion. Each successive constraint involves the explicit formula for the higher derivatives of 
coupled with the extra constraint that 

|@(n)|. = { Q (n+1)}. 

We conjecture that each constraint can be expressed as a commutator equation in terms of , 
H and the derivatives of H, in analogy to the formulas that we have found for the first three 
constraints. This project will continue with a deeper algebraic study of the constraints and their 
physical meanings. 



8 Appendix - Einstein's Equations and the Bianchi Identity 

The purpose of this section is to show how the Bianchi identity (see below for its definition) 
appears in the context of non-commutative worlds. The Bianchi identity is a crucial mathematical 
ingredient in general relativity. We shall begin with a quick review of the mathematical structure 
of general relativity (see for example [6]) and then turn to the context of non-commutative worlds. 
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The basic tensor in Einstein's theory of general relativity is 

G ab = R ab _ \ R ab 

2 y 

where R ab is the Ricci tensor and R the scalar curvature. The Ricci tensor and the scalar curvature 
are both obtained by contraction from the Riemann curvature tensor R bcd with R a b = R c ahc , R ab = 
g at g b i Rij , and R = g l ^Rij . Because the Einstein tensor G ab has vanishing divergence, it is a prime 
candidate to be proportional to the energy momentum tensor T^ u . The Einstein field equations 
are 

2 y 

The reader may wish to recall that the Riemann tensor is obtained from the commutator of a 
covariant derivative associated with the Levi-Civita connection F l jk = (F k )'j (built from the 
space-time metric gij). One has 

or 

\. b = v h A = 3 b \ - T b \ 

for a vector field A. With 

Rio = [Vi, V,] = djTi - diYj + [Fi, Fj], 

one has 

R-bcd = \Rcd)b- 

(Here R C( i is not the Ricci tensor. It is the Riemann tensor with two internal indices hidden from 
sight.) 

One way to understand the mathematical source of the Einstein tensor, and the vanishing of 
its divergence, is to see it as a contraction of the Bianchi identity for the Riemann tensor. The 
Bianchi identity states 

R-bcd-.e + R-bde-.c + ^bec-.d = 

where the index after the colon indicates the covariant derivative. Note also that this can be 
written in the form 

(Rcd:e)b + {Rde-.cjb + (Rec:d)b = 0- 

The Bianchi identity is a consequence of local properties of the Levi-Civita connection and con- 
sequent symmetries of the Riemann tensor. One relevant symmetry of the Riemann tensor is the 
equation Rfa = -R% dc . 

We will not give a classical derivation of the Bianchi identity here, but it is instructive to see 
how its contraction leads to the Einstein tensor. To this end, note that we can contract the Bianchi 
identity to 

^bca-.e ' ^bae-.c ' - rt 6ec:o u 
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which, in the light of the above definition of the Ricci tensor and the symmetries of the Riemann 
tensor is the same as 

Rbc-.e ~ Rbe-.c + Rbec.a = 0- 

Contract this tensor equation once more to obtain 

Rbc.b — Rbb-.c + Rbbc.a = 0) 

and raise indices 

Rc:b ~ R-c + Rbc-.a = 0- 

Further symmetry gives 

n bc:a — n cb:a ~ n c:a ~ n c:b- 

Hence we have 

2R c: b ~ R-c = 0, 

which is equivalent to the equation 

(R b c ~ 2^c):6 = G b c .b = 0. 

From this we conclude that G b £ = 0. The Einstein tensor has appeared on the stage with vanishing 
divergence, courtesy of the Bianchi identity! 

Bianchi Identity and Jacobi Identity. Now lets turn to the context of non-commutative worlds. 
We have infinitely many possible convariant derivatives, all of the form 

F-. a = V a F= [F,N a ] 

for some iV a elements in the non-commutative world. Choose any such covariant derivative. 
Then, as in the introduction to this paper, we have the curvature 

Ri^lNuNj] 

that represents the commutator of the covariant derivative with itself in the sense that [Vj, Vj]F = 
[[iVj, Nj] , F] . Note that Rij is not a Ricci tensor, but rather the indication of the external structure 
of the curvature without any particular choice of linear representation (as is given in the classical 
case as described above). We then have the Jacobi identity 

[[N a , N b },N c ] + [[JV C , N a ],N b ] + [[N b , Ay, N a ] = 0. 
Writing the Jacobi identity in terms of curvature and covariant differention we have 

Rab.c ~\~ Rca.b Rbc:a- 

Thus in a non-commutative world, every covariant derivative satisfies its own Bianchi identity. 
This gives an impetus to study general relativity in non-commutative worlds by looking for co- 
variant derivatives that satisfy the symmetries of the Riemann tensor and link with a metric in an 
appropriate way. We have only begun this aspect of the investigation. The point of this section 
has been to show the intimate relationship between the Bianchi idenity and the Jacobi identity 
that is revealed in the context of non-commutative worlds. 
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9 Philosophical Appendix 



The purpose of this appendix is to point to a way of thinking about the relationship of mathemat- 
ics, physics, persons, and observations that underlies the approach taken in this paper. We be- 
gan constructions motivating non-commutativity by considering sequences of actions • • ■ DCBA 
written from right to left so that they could be applied to an actant X in the order ■ ■ ■ DCBAX = 
■ ■ ■ (D(C(B(AX))) ■ ■ ■ . The sequence of events A, B,C, D, • ■ ■ was conceptualized as a tempo- 
ral order, with the events themselves happening at levels or frames of successive "space". There 
is no ambient coordinate space, nor is there any continuum of time. All that is given is the possi- 
bility of structure at any given moment, and the possibility of distinguishing structures from one 
moment to the next. In this light the formula DX = [X, J] = XJ — JX = J(X' — X) connotes 
a symbolic representation of the measurment of a difference across one time interval, nothing 
more. In other words DX represents a difference taken across a background difference (the time 
step). Once the pandora's box of measuring such differences has been opened, we are subject to 
the multiplicities of forms of difference V^X = [X, K] , their non-commutativity among them- 
selves, the notion of a flat background that has the formal appearance of quantum mechanics, 
the emergence of abstract curvature and formal gauge fields. All this occurs in these calculi of 
differences prior to the emergence of differential geometry or topology or even the notion of 
linear superposition of states (so important to quantum mechanics). Note that in this algebraic 
patterning each algebra element X is an actant (can be acted upon) and an actor (via the operator 
Vx). In Lie algebras, this is the relationship between the algebra and its adjoint representation 
that makes each element of the algebra into a representor for that algebra by exactly the formula 
adjA(X) = [A, X] = — Va{X) that we have identified as a formal difference or derivative, a 
generator for a calculus of differences. 

The precursor and conceptual background of our particular formalism is therefore the con- 
cept of discrimination, the idea of a distinction. A key work in relation to that concept is the 
book "Laws of Form" by G. Spencer-Brown G8l in which is set out a calculus of distinction of 
maximal simplicity and generality. In that calculus a mark (denoted here by a bracket < >) rep- 
resents a distinction and is seen to be a distinction between inside and outside. In this elemental 
mathematics there is no distinction except the one that we draw between the mathematician and 
the operator in the formal system as sign/symbol/interpretant. This gives full responsibility to the 
mathematician to draw the boundaries between the formal system as physical interaction and the 
formal system as symbolic entity and the formal system as Platonic conceptual form. In making a 
mathematics of distinction, the mathematician tells a story to himself/herself about the creation of 
a world. Spencer-Brown's iconic mathematics can be extended to contact any mathematics, and 
when this happens that mathematics is transformed into a personal creation of the mathematician 
who uses it. In a similar (but to a mathematician) darker way, the physicist is intimately bound to 
the physical reality that he studies. 

We could have begun this paper with the the Spencer-Brown mark as bracket: < > . This 
empty bracket is seen to make a distinction between inside and outside. In order for that to 
occur the bracket has to become a process in the perception of someone. It has to leave whatever 
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objective existence or potentiality it has alone (all one) and become the locus or nexus of an idea 
in a perceiving mind. As such it is stabilized by that perception/creation and becomes really a 
solution to{< >}=< > where the curly bracket (the form of perception) is in the first place 
identical to the mark < >, and then distinguished from it by the act of distinguishing world and 
perceiver. It is within this cleft of the infinite recursive and the finite 

< >={< >} = {{< >}} = {{{< >}}} = •■■ = {{{{{{•■■}}}}}} 

that the objectivity of mathematics/physics (they are not different in the cleft) arises. All the 
rest of mathematics or calculus of brackets needs come forth for the observer in the same way. 
Through that interaction there is the possibility of a deep dialogue of many levels, a dialogue 
where it is seen that mathematics and physics develop in parallel, each describing the same 
boundary from opposite sides. That boundary is the imaginary boundary between the inner and 
outer worlds. 
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